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1 Introduction 

Investigating the decoupling limit of massive gravity-related theories (see [1, 2] for reviews) 
has turned out to be a very fruitful exercise over the past few years. The decoupling limit, 
by which we mean the scaling limit of these theories isolating non-linear interactions between 
physical degrees of freedom ( dof’s ) suppressed by the smallest energy scale, reveals these 
dof's and their low-energy interactions as well as strong coupling scales for the theories in 
question. Following in the footsteps of work investigating the decoupling limit of ghost-free 
(dRGT) Massive Gravity (see [3-5] for the model itself and [3, 6-9] and references therein 
for its decoupling limit interactions) and (Hassan-Rosen) Bigravity (see [10-12] for the model 
itself and [13] for investigations of its decoupling limit), we here probe the decoupling limit of 
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(Hinterbichler-Rosen) Multi-Gravity theories [14] for the first time. For constraint analyses 
relevant for these models see [14, 15] and for related multi-gravity work also see [16-20]. 1 

One of the ways in which investigations of the decoupling limit have been very fruitful 
is by uncovering a set of dualities non-trivially relating (seemingly different, but ultimately 
physically equivalent) scalar [13, 22, 23] and vector [23] field theories via non-local, invertible 
field re-definitions. In particular these dualities can sometimes be useful in relating strongly to 
weakly coupled models. Scalar field dualities have been uncovered in this way by investigating 
helicity-0 mode interactions in the decoupling limit of Bigravity [13] . These are described by 
Galileon interactions [24] (also see [25, 26]) for the helicity-0 scalar dof (as is also the case 
for massive gravity) and the dualities relating these interactions can also be understood as 
abstracted duality transformations [ 22 ] resulting from transformations on the coset space of 
GAL(d, l)/SO(d— 1,1) [27]. The duality can also be extended to ‘generalized galileons’ [23]. 

The existence of these dualities is deeply linked to the non-uniqueness of the way in which 
gauge dof's can be added to the theory in order to restore full diffeonrorphism invariance - 
a feature that is made explicit in the link held formulation for Multi-Gravity theories [28— 
32]. Using this formulation and the decoupling limit interactions for classes of Multi-Gravity 
models, in this paper we extend the (single) held Galileon dualities hitherto discovered to 
Multi-Galileon dualities and show how they can be understood both as inherited from diffeo- 
morphism invariance and as abstracted duality transformations. 

Outline: This paper is organised as follows. In section 3 we review decoupling limit in¬ 
teractions for massive and bigravity and extend this to a number of ghost-free multi-gravity 
setups. In section 4 we then discuss Galileon dualities in two ways: Firstly as a direct con¬ 
sequence of (diffeomorphism) gauge invariance in ghost-free Bigravity models and secondly 
as an abstract duality transformation. We then use the insight gained in this way to derive 
a set of multi-galileon dualities in section 5, where we also show how this can be used to 
demonstrate the healthiness of interactions between helicity-0 in multi-gravity models. Fi¬ 
nally we discuss how matter can be coupled in multi-gravity theories and how this relates to 
the duality picture in section 6 , before concluding in section 7. 

Conventions: Throughout this paper we use the following conventions. D refers to the 
number of spacetime dimensions and we use Greek letters /r, is, .. . as well as lower case Latin 
letters a, b, ... to denote spacetime indices, which are raised and lowered with the Minkowski 
metric r]^ u unless explicitly stated otherwise. We write the completely anti-symmetric epsilon 
symbol as e and define it such that eoi 2 ---D = 1 regardless of the signature of the metric or 
the position (up/down) of indices (hence e 012 " D = 6012 - D = !)• 

1 Also, for a discussion of some of the as yet unresolved issues in the general field of massive gravity(-ies) 
and whether they are problematic for the field see e.g. [21] and the counterarguments in [2], 
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2 Degrees of freedom, equations of motion and the absence of ghosts 
2.1 Degrees of freedom and diffeomorphism invariance 

General relativity, a theory of a massless spin-2 particle, propagates two dof’s (around flat 
space and in the absence of a coupling to matter). A massive spin-2 field on the other hand 
propagates 5 dof’s. When building a massive gravity action, however, it is typically not ob¬ 
vious exactly how many dof’s propagate and what their dynamics is. The Stiickelberg trick 
is useful in this circumstance, mapping a given action into another dynamically equivalent 
action, which contains more fields and more symmetry. In essence the Stiickelberg trick is 
therefore the inverse of gauge-fixing an action, taking an action without a given symmetry 
and restoring it via the addition of gauge dof’s. For a more comprehensive review see [32] - 
here we will be content to only give the bare recipe for performing the Stiickelberg procedure. 2 

Massive Gravity. Consider a generic local and Lorentz-invariant theory of massive 
gravity formulated in terms of a ‘metric’ g and a non-dynamical, flat background metric r/ 

S\ = M°~ 2 j d D x^fjR + m 2 Mp l ~ 2 J d D x^~gV (g~ x r{) . (2.1) 

The potential V (<? -1 ^) breaks diffeomorphism invariance, which is just to say that it leads 
to the propagation of additional dof’s. We may restore full diffeomorphism invariance by 
performing either of the two following Stiickelberg replacements patterned after the symmetry 
we are restoring and at the expense of introducing extra (gauge) fields 

9»u d^Y a d v Y^g aP [F(x)], (2.2) 

ri, UJ -»• c^WUTV/3- (2-3) 

Y a and Y° are the Stiickelberg fields encoding the extra gauge field content. That we have 

a choice of how to introduce these fields is a simple consequence of the fact that there is no 
unique procedure to introduce the redundancy associated with gauge fields when adding these 
fields to the action in the process of restoring diffeomorphism symmetry. Since in massive 
gravity we have a non-dynamical metric rj, the result of using (2.2) is an object which is 
invariant under the restored single copy of diffeomorphism invariance, rather than a tensor 
and so is more complicated to analyse than the result of using (2.3) - for details see [1]. We 
can now expand the Stiickelberg fields via 

Y 11 -> x 11 + + &*<!>, 

-> x^ + + &*ir, (2.4) 

respectively. Y^ = x^ 1 and Y 11 = x IJ/ correspond to fixing the gauge fields to ‘unitary gauge’. 
In expanding the Stiickelberg fields we have also introduced an additional 17(1) symmetry 

2 A1so see [33] for recent work on ‘covariant Stiickelberg analyses’. 
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(effectively performing a second Stiickelberg trick) with the associated fields n and </> re¬ 
spectively. In the A 3 decoupling limit A/B will describe the helicity-1 modes and 7 r/cj) the 
helicity-0 modes of the spin-2 fields in question. An important point is that the replacement 
(2.3) is finite order in A and 7 r, whereas (2.2) is not due to the dependence on g a p [y(x)] (a 
straightforward way to see this is to Taylor-expand g a p [T(x)] ~ see [1] )- This makes identify¬ 
ing the interactions between different helicity modes very straightforward when (2.3) is used 
to restore diffeomorphism invariance in the action . 3 

Bi- and Multi-Gravity. Now consider a bi- or general multi-gravity model of the 
following type 

N ' /• _ p 

S 2 = ^2 M pf 2 / d D x^/-g { i)R [g {i) ] + m 2 M^ 2 / d D x^J-g {l) V (c/ (1) ,..., fl'(jv))) , (2.5) 

i =1 •' 

where we have N dynamical spin-2 fields described by gjll , where bracketed indices ( i) are 
label indices only summed over if there is an explicit sum in what follows (i.e. there is no 
Einstein convention for label indices throughout this paper). The potential is a Lorentz-scalar 
built from the different g^ and their inverses and all N fields are dynamical in contrast to the 
massive gravity model above, which was really a bigravity theory with one non-dynamical 
(fixed) field 77 . 4 Note that the choice of metric determinant in front of the potential can 
be absorbed into the definition of the potential and is therefore arbitrary. Let us focus on 
the bigravity version of (2.5). The potential V now generically breaks the two copies of 
diffeomorphism invariance present at the level of kinetic interactions for g^\ and g^) down to 
the diagonal subgroup. The broken copy can again be restored via either of the two following 
replacements 



where Y and Y can be expanded as before. Note that we no longer have the luxury of being 
able to introduce the Stiickelberg fields via a non-dynamical field rj as before for theories 
of the type (2.5), so these replacements are infinite order in the Stiickelberg fields and in 
derivatives. As we shall see the fact that we can choose between two such replacements and 
hence two equivalent ways of capturing the helicity-0 mode ( 7 r and </>) here is intrinsically 
linked to the existence of Galileon dualities. 

2.2 Higher-derivative equations of motion 

Having performed the Stiickelberg replacement discussed above one can read off the interac¬ 
tions of the different helicity modes and their interaction scales. The first check one would 

3 Note that there is of course also a choice as to what symmetry one would like to restore, if any - see e.g. 
[34] where instead of full diffeomorphisms, transverse diffeomorphisms and a Weyl symmetry are restored. 

4 We have also assumed that the Planck masses for all spin-2 fields are the same and that the kinetic sector 
is solely made up from Einstein-Hilbert kinetic terms - for details on more complex kinetic sector interactions 
see [35-40]. 
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(a) Bigravity 


(b) (Trimetric) line theory 


(c) Star theory 


Figure 1. Theory graph representations of different types of Multi-Gravity theories. Each node 
represents a metric gm and a double-circled node denotes the presents of an Einstein-Hilbert kinetic 
interaction yJ—g(i)R[g(i)\- Lines connecting two nodes are bigravity-like interaction terms between the 
two fields associated with those nodes (note that we do not include interactions between more than 
two metrics in this paper). Shaded nodes are coupled to matter (minimally in the cases shown here). 
See [40] for further details. The graphs shown depict: a) Hassan-Rosen Bigravity, b) A particular 
trimetric ‘line theory’ as considered in sections 3.3 and 5, with no coupling to matter specified, c) A 
‘star theory’ as considered in section 6, where several outer nodes are connected to a single central 
node, but not to each other. In the particular case shown here, only the central node is (minimally) 
coupled to matter. 


typically like to carry out, is to see whether the interactions found for the model under con¬ 
sideration are consistent - in particular: whether they are free of ghost-like instabilities. With 
an eye on the multi-gravity models we will consider in more depth later on, let us focus on the 
interaction between two scalars n and (f> at some scale A (the scalars mimic helicity-0 modes 
and A mimics the decoupling limit scale). When can we say these interactions are healthy 
(at least up to the scale A)? 

A multi-galileon. Let us start with a familiar Lagrangian, namely that of a particular 
cubic multi-galileon (cf. [41-43]). 

Ci = ^7rD7r + ^</>□</> + ^7r (d b d a <j)d b d a (j) - (D^) 2 ) . (2.7) 

The equations of motion ( eoms ) for this Lagrangian are 

£</, = □</>+ {dbd a 'nd b d a (t) - n^jDnj = 0 , 

£* = j (d b d a (/)d b d a (t> - (n^ 2 ) = o. (2.8) 

These eoms are manifestly second order and the linear kinetic term has the correct sign. This 
is a quick diagnostic to show that the system is free of (Ostrogradsky) ghost-like instabilities. 
If the eoms have dependence on higher-order derivatives, in the absence of additional con¬ 
straints/degeneracies this will lead to propagation of an extra ghostly dof commonly called 







an Ostrogradsky ghost [44, 45]. 


A sick theory. A quick example of such a ghostly, higher order in derivatives Lagrangian 
is given by a small modification of £], where we have broken the antisymmetric structure of 
the cubic interaction term 

£2 = -vrDvr + -</>□</> + — 7rn7rD(/). (2-9) 

The corns now are 

4 = U<t> + j ((Dtt) 2 + 2d a 7r<9 a D7r + 7rD 2 7r) = 0, 

£ n = D7r + (2D0D7T + 2d a ird a O(j) + 7rD 2 0) = 0, (2-10) 

and they are explicitly higher-order. In particular, the first eom depends on up to fourth 
derivatives of n, whereas the second eom depends on up to fourth derivatives of qb. Con¬ 
sequently more than two initial conditions need to be specified per held and the extra dof 
encoded in these additional initial conditions is an Ostrogradsky ghost. Note that, if an in¬ 
vertible, non-linear held re-dehnition can be found, which removes interactions at this order 
in the helds, the addition of higher order (in the helds) interactions to the Lagrangian can 
prevent the appearance of a ghost at the scale A. 


Healthy higher-derivative eoms. Finally let us consider a case, where the eoms appear 
higher-order, yet are ghost-free due to the existence of an extra constraint. The following is 
a straightforward example of this situation, which we have lifted from [46] 

£3 = ^7rD7r + ^ </>□</> + ^n-nd a d b (f)d a d b (j) + ^^d a d b (l)d a d b (j)n(dcd d (j)d c d d (j)) ( 2 . 11 ) 

The equations of motion here are 


St = n<i> + ^d a d b (d a d b cp □*- + (<w<9 d c>v) 


= 0, 


£tt = Dtt + J^\3(d d d e (j)d d d e (j)) = 0, 


( 2 . 12 ) 


which appear to be higher order in derivatives. However, note the explicit dependence of the 
first eom on the second, which is the way in which the extra constraint manifests itself here. 
By writing £^ = £ r j ) — 2/A 5 d a d b (d a d b 4>£ n ) and £ n = £ n , we can instead write the eoms for cj) 
and 7 r as 


4 W 1 ] =□</> = 0 , 

£tt [tt 11 , (t > H , </> m , 4> n ] = Dvr + ^U{d d d e (j)d d d e (j)) = 0, (2.13) 

where we have made the dependence on derivatives explicit, i.e. £ n [k 11 ,...] denotes that 
this eom depends on second derivatives acting on 7r and so on. It is then clear that cj) can 
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be solved for in terms of just two initial conditions. £- n is higher-derivative in nature, but 
the higher-derivative dependence is restricted to </>, for which we have already solved. So 
7 T can also be solved for in terms of two initial conditions and the whole system is secretly 
second-order with a well-defined Cauchy problem. Schematically we have derived a final set 
of eoms, £^,£( 1 ), where = £<p and £ n was obtained by substituting the solution of into £ n . 
Consequently we have two final eoms £$ [4> n ] ,£ n [vr 7/ , 4> n ]. The same conclusion could have 
been reached by performing an explicit constraint analysis for this model in the Hamiltonian 
picture. In general, if there is any linear combination of the eoms of a system and derivatives 
of the eoms that explicitly allows all the fields to be solved in terms of two initial conditions for 
each field (here we are assuming we are dealing with scalars, obviously this count is modified 
for other tensors), then we have a second-order system free from Ostrogradsky ghosts. Note 
that the models proposed in [47] as well as the recently proposed ‘beyond Horndeski’ theories 
[48, 49] are precisely of this type as well (albeit being theories of a scalar and a metric tensor, 
instead of multiple scalars as discussed here, in the case of ‘beyond Horndeski’ theories). 


2.3 Some notation 

Let us quickly summarise some conventions used throughout the remainder of this paper. We 
will reserve the letters it and <for Stiickelberg scalars. In order to keep notation concise, we 
will typically simply denote partial derivatives acting on these scalars by indices, i.e. 

<:; b a n = d ai ...d an d b i...d b "7r (2.i4) 


Both Latin and Greek letters denote space-time indices unless stated otherwise. Bracketed 
indices as in 7 Tm are label indices, e.g. labelling different helicity -0 modes in multi-gravity 
setups. Label indices are never summed over. In order to connect with notation used in the 
literature, we will sometimes denote 7 r 7 by n^, 

H'i(x) = ^V^dadMx) , (2.15) 

i.e. we will reserve n to denote this particular two derivative function of n, especially when 
using index-free notation - see below. A is a mass scale, which in the 4D massive gravity 
context will turn out to be A 3 = ( m 2 Mpi ) 1//3 . For a scalar 7 r, at each order n, there is a 
unique total derivative combination given by 


£ ™( n >= 


A 

ai 


• 7TCT 


(2.16) 


where we have defined a tensor separately anti-symmetric in its indices a\... a n and 

Pi ... f3 n in terms of the totally antisymmetric tensor e via 


rOi ...a n _ 1 

Ifr-M ~ (. D-n)\ 




£fll-../3nAl...A£)_n 


(2.17) 


Equivalently to £™(n), and in order to connect with the notation used by [13, 22, 23], we 
may also define the characteristic polynomials C// n ) of a matrix AI in the following index-free 



way 


U(n){M) = eeM n r 1 D - n , (2.18) 

where f/( n ) (II) = (D — n)!£™(II). To give an explicit example, for [/( 2 ) in 4D this means 

17 (2) (n) = e^e^U^U^pcri^ 

— 9U" 1 “ 2 “ 3a4 TT^ 1 TT^ 2 A'® 3 A^ 4 

iia 2 °«3 0 O4 

= 2£™(II). (2.19) 


In addition we can define the tensor Tf^(II), which is the unique symmetric and identically 
conserved function of II at a given order n 


%]'.(n) = «f“i n 2 ■ ■ • n h < 22 °) 

We will also find it useful to define the related tensors X^ v (Jl) and X^ v {Ii), which mix all 
orders in n 


X^(U) 

X^(U) 


l ^ 1 _ l 

2 ^ (D - 1 


Pn 


2 ' (D — 1 — n)!n! 

n=o v ' 


e »-e v -{rt + Tl) n r 1 


n^D—l—n 


i D 

E 


A, 


D— m^m —1 


2 (£) _ rn)\(m — 1)! 

m= 1 


e^-e^-^ + nr"^ 


( 2 . 21 ) 


where the /3 n are constant coefficients. Note that X and X can straightforwardly be mapped 
into each another with the replacement m = D — n - they are essentially the same object 
with remapped coefficients 5 . Characteristic polynomials , the ‘X-tensors’ defined above and 
total derivative combinations at a given order are related by 


V,uX^(U) 

n) 

V^(U) 

yT)( n ) 


(D-l-*)!^ 0 ^’ 

1 Ar 


2 ^ (D — 1 — n)\n\ 
n= o v J 

1 Pm 

I 

1 5 


2 ' (D — m)\(m — 

171 =0 X 7 V 


U(n){ji + n), 

Yyi^M^ + n), 


£™(n). 


n + 1 hn Mt/ 


( 2 . 22 ) 


3 The decoupling limit of Massive, Bi-, and Multi-gravity 

In this section we will review the decoupling limit interactions for massive and bigravity and 
derive analogous results for a class of multi-gravity models. 

5 Note that our X is essentially the same object as Y defined by [13]. We choose our notation in order to 
emphasise that X and X denote identical interactions up to constant coefficients and since this definition will 
be rather useful in the general multi-gravity treatment of section 5. 
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3.1 Massive gravity 

The ghost-free dRGT massive gravity model [3-5] can be written down as 


S = d D x 


D 


M»- 2 ^gR[g\+m 2 M, 


D—2 
PI 


V~9 X @ n u ^(Vg lr i) 


n =0 


(3.1) 


where the j3 n are constant coefficients, U n denotes an elementary symmetric polynomial of 
order n and (in index-free matrix notation) y/g~hj = M, where M is any real matrix satisfy¬ 
ing M 2 = We now perform the Stiickelberg procedure, restoring full diffeomorphism 

invariance as shown above (we here choose (2.3)) and making interactions of the different 
helicity modes explicit in the process. Expanding the metric around a flat background, 
Qfiu = Vfiu + h^, one then finds that interactions involving n n fields 7 r, ua fields A^ and rih 
fields h^, carry an interaction scale 


Aa = 



_ \ 1/A 

m A_1 I = m 



A 


3n n + 2 ua + nh - 4 
n-K + ua + n h - 2 


(3.2) 


The least suppressed interactions that survive in the dRGT model carry a scale A 3 .' We can 
isolate interactions at this scale by taking the A 3 decoupling limit 6 * 8 


m 


0 , 


Mpi —> 00, 


A 3 fixed. (3-3) 


The helicity-1 mode A does not appear linearly in the action, so we can consistently set it to 
zero, which we will do for the remainder of this paper . 9 

In the A 3 decoupling limit the helicity-2 and helicity-0 modes are described by h and 7 r 
respectively - they are now governed by the following action 


S 



-4 + Y X <*(n)h'“'X$(n) 

n =1 


(3.4) 


6 Note that the existence of such a real square root is a non-trivial requirement, related to the symmetric 
vielbein condition [15]. 

'In 4D this is a mass scale. 

8 Note that in this limit the dynamics of the helicity-0 mode is captured by the Stiickelberg scalar n, the 
helicity-1 modes are captured by the A M and the helicity-2 modes are captured by whereas beyond the 
decoupling limit (i.e. for more suppressed interaction terms), this is not the case [50]. This is beyond the 
scope of this paper, however, and we restrict ourselves to a discussion of the decoupling limit phenomenology 
throughout. 

9 In fact, as long as matter is minimally coupled to the metric, A never gets sourced. For details on helicity-1 
interactions in massive gravity see [7, 8]. 
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where the are constant coefficients and we have ignored the coupling to matter. 10 
can now demix the interactions between h and ir via the following field redefinitions 

2ctn ) Ag 

hjj,v ^ hfti, H — —TtTliivi 
hfti/ > — 2ct(2)AgTr^TTj,. 


We 


(3.5) 

(3.6) 


The first one, a linearised conformal transformation, demixes the different helicity modes 
at linear order, eliminating the h! iv X^,) (7r) mixing. Subsequently applying the second field 
re-definition also eliminates the cubic h ,iu XfJ (-/r) mixing - demixing at higher orders cannot 
be achieved with a local field re-definition. After these replacements, the helicity-0 self¬ 
interactions at cubic, quartic and quintic order (in 4D these are the only non-vanishing terms) 


are given by 

5 ( 3 ) ~ J d D X TT a TT a TT^, (3.7) 

5(4) ~ J d D x ir a ir a (n%TT° - , (3.8) 

5(5) ~ J d D x n a ir a (3ir d d ir bc ir bc - 2ir d c -K db ir bc - ir d d ir b b ir c c ^J . (3.9) 

Note that S 3 and 54 are already present after the first field redefinition (3.5), whereas 5s 
only appears after (3.6) has been performed. 11 The corresponding equations of motion are 
5(j) = 0, where 

£(3) = KA - = (X) * 2)! C/(2) ( n )’ ( 3 ‘ 10 ) 

5(4) = 3vr a a 7T bc 7r bc - 2vr a c 7r ab 7r bc - Tr a a Tr b b ir c c = — -j-^ [/^(II), (3.11) 

5( 5 ) = 8TT a a TT b d TT bC TT cd - §’K a C 'K ab TT b d 'K cd + 8TT ab TT ab TT cd TT Cd (3.12) 

- 6TT a a TT b b TT cd 7T Cd + TT a a TT b b TT C c 7r d d = — — — f7( 4 )(II). 


As is well-known, these terms are precisely the Galileon interactions [24] and their corre¬ 
sponding eoms (where we avoid any integration-by-parts ambiguities). Consequently we may 
write the action for the helicity-0 mode in the decoupling limit as 

S = f d D xJ2c(i)irU n (Il). (3.13) 

i 

Here all higher order equations of motion for ir vanish and we can in fact look at each order 
in the eoms separately to confirm that they are (purely) 2nd order in derivatives of 7 r at all 
orders. Consequently 7r does not give rise to any ghost-like (Ostrogradsky) dof. 

lll Note that -X^t) ( 77 ) and higher orders in X vanish identically in 4D. So the hn couplings shown in (3.4) are 
all the non-vanishing h^ 1 ' xffi (it) terms in 4D. 

11 This is also a quick way of seeing that there is only a two-parameter family of solutions present in the 
decoupling limit - a statement which is also true for the full theory. 
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3.2 Bigravity 

The Hassan-Rosen bigravity model [10-12], essentially the extension of the dRGT massive 
gravity model in which both metrics are dynamical, is described by the following action 


S = 


d D x 


D 


Mj 


D—2 


PI 


V^9 R[g\ +m£" 2 v/=/ R[f] +m 2 M°- 2 ^r g £ 3 n U n {y/g=^f) 


n =0 


, (3.14) 


where we have assumed that the Planck masses for both metrics are identical - a generalisation 
is straightforward. As before we perform the Stiickelberg replacement, ignore the helicity-1 
modes for the time being and focus on the interactions of the helicity -2 and helicity -0 modes 
in the decoupling limit . 12 Due to a symmetry of the elementary symmetric polynomials 
namely ^gU (n) ( \fg~ Y f) = f~ l g), these can be written [13] 


5 = 


d s 5 


~4 (i)' » v 4 


(3.15) 


where cr is a field ‘dual’ to it and can be implicitly defined via 13 


x + da = (x + Ott) \ 


(3.17) 


and in analogy to our definition of II we have also defined £ via Yf v (x) = d fl d u a(x) /A|. We 
can then iteratively solve for a in terms of it and find (to quintic order in n) 


_ _ | ^ _ ,_a ^ i ^ ^ 

CJ — 7T i - q 7T/ti7T r, 7T 7T TTnh ~i . n 7T 7T TT/i TT/)^ — i 7T 7T 7T TTnhc 

2 A| 2A| 2A§ 6 A| 

^ c_ d _ ^ ^abc„ d_ ^ _a_6^c_d_ , 

2 A12 71 ^ Ti 2 A 12 ^ ^ ^ Ti 24A 12 ^ ^ ^ ^ Ttabcd T C/(7T ). 


(3.18) 


We will look into the relation between cr and 7 r in more detail in section 4 (in particular a 
full solution is given in (4.13)), but looking ahead it is worth pointing out that cr is related 
to 7 r via a non-local, but invertible, field redefinition [13, 22], This is important because such 
a field redefinition preserves the number of dof’s (in this case the dof described by n). 


Diagonalising kinetic terms. We now demix the action at linear order (for full details 
see [32]). There are a number of possible field re-definitions we may choose to do so, two of 
which are 

Mi : h$ -> h$ + civrr/, h$ -> h$ + c 2 m 7 , 

M 2 : h$ -> h$ + civrr/, hg) -> hg) + c 2 tu]. (3.19) 

12 Since the two Planck masses are the same, this limit is the same as before, otherwise different decoupling 
limits are possible [32]. 

13 We can also define cr via an implicit function Z at this point [13] 

Z a (x b + gg d b n(x)) = x a , Z a (x) = x a + ~^d a a(x). (3.16) 

■^3 A3 

We will discuss the origin and form of a in more detail in the next section. 
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As far as demixing at linear order goes, the difference between Mi and M 2 is irrelevant, since 
to first order a{ it) = —n, but we will see that .Mi makes the number of propagating dof’s 
and their interactions more explicit. When applying Mi we obtain the following schematic 
action for the pure helicity-0 piece post-demixing (and still in the decoupling limit) 

Si = J d D x Y, {c(i)KU n ( n) + d (i) <rtf n (E)) , (3.20) 

i 

i.e. we have an explicitly second-order galileon-like term for 7 r and one for a. However, a 
and 7T describe the same dof's and we should therefore really express the action in terms of 
just one of the two. What is not trivial to see is that the this action also only propagates a 
single scalar dof when a is expressed in terms of 7r or vice versa. Naively one may expect that 
the then higher derivative nature of the action would lead to the appearance of additional 
propagating dof’s. 11 When applying M 2 , on the other hand, we obtain 

S 2 = f d D xY (c W 7r^ n (n) + d W 7rC4(£)) , (3.21) 

i 

where it is now no longer clear that the second term describes a Galileon-like interaction and 
the true number of propagating dof’s is not obvious. 

Equations of motion. We now again move on to examine the eom for 7r, where we 
directly express a in terms of ir via (3.18). We already know that the 7rt/ n (n) piece in (3.20) 
will contribute second-order Galileon terms to the eom. In addition from the aU n (Yi) in (3.20) 
we obtain the following contributions to the eom when varying with respect to 7r 

£( 1 ) = 0 , 

8(2) = Pi* a a , 

8(3) = Pli^a^b - TTabTT ab ) + P 2 ( §7r a6 7r“ 6 - ^TT a a TT b b ), 

8(4) = (dl{TT a C TT ab TTbc - |vr a a 7T fec 7r feC + 5 Tr^TT^c) 

+ a^bcTT bC - \n a bK C C ~ 37T a C TT ab TT bc ) 

+ /3 3 (27r a C 7T a6 7 T bc - 3-K a a 7T bc TT bc + !T a a lT b b lT C c ), (3.22) 

where 8 (^ denotes the contribution to the eom coming from the i-th order (in n) piece in the 
action. 8 ^ is consequently (i — l)-th order in tt. We have also omitted quintic order (this 
can be found in appendix A) and set A 3 = 1 to avoid clutter. We can immediately see that 
all the contributions are second order in derivatives and this remains true at quintic order, 
i.e. the Galilean symmetry is manifest. Importantly, up to quintic order in the fields at least, 
there is therefore no Ostrogradsky instability arising from helicity-0 interactions in ghost-free 
bigravity. 

14 In this subsection we will directly inspect the eoms and deduce the number of propagating dof’s that way. 
Independently the bigravity duality picture as discovered by [13, 22] and discussed in section 4 will also show 
independently that no additional dof propagates due to the apparent higher-order nature of the action. 
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Suppose we had demixed differently with A4 2 and consequently obtained (3.21) instead. 
The extra (not obviously of Galileon form) contribution to the eom now comes from the 
irU n (Ti) piece in (3.21). This results in the following contribution when varying with respect 
to 7 r 

£(i) = o, 

£ ( 2) = Avr a a , 

£(3) = - l^ab^ ab ) + P 2 {\n ab 'K ab ~ f 7T a a 7r\), 

c _ o(l— c_a&_ 1 a _ be , 1 a b c 

£(4) — Pl\ 2 ^ a ^ ^be 4 ^ a^bc^ “ 1 “ 4 ^ 6 7T c 

+ TT a ir bc ir abc + ^ir a TT b b 7T a C c + \-K a TT a b ^b C c + 57 r a 7 T b 7 r afe c c ) 

+ /3 2 (37r a a 7r fec 7T 6c - 27r a c 7r a6 7r fec - Tr a a ir b b n c c ) 

+ /3 3 (2ir a c Tr ab ir bc - 37 r a a 7 T fec 7 r 6c + n a a Tr b b ir c c ). (3.23) 

Up to cubic order in 7 r and ignoring irrelevant overall coefficients these contributions to the 
eom are identical to those derived from (3.20). However, at quartic order they become different 
and noticeably we now have a higher-derivative dependence (the second line of £( 4 ) contains 
all these terms). In particular these terms all only depend on j3i. If we continue to quintic 
order (see appendix A) we find that there is higher-derivative dependence via f3 2 now as well. 
This is a direct consequence of the fact that the two different demixing procedures agree up to 
linear order, but diverge at higher orders in the field. The higher derivative dependence of the 
second set of contributions (3.23) is inherited from the action’s dependence on a, cf. (3.18), 
which is non-locally related to 7 r, i.e. infinite order in derivatives. The fact that the action now 
non-locally depends on 7 r means that Ostrogradsky counting is not(!) appropriate for the full 
theory - the action is degenerate since one cannot isolate a highest-order-in-derivatives term 
[44, 45]. The fact that the first set of contributions (3.22), obtained via a different demixing 
procedure, is purely second-order is already a direct hint that a non-local, invertible field 
re-definition exists for 7 r in the second set of contribution (3.23), which cures the dependence 
on higher derivatives. (3.23) is therefore another example of an eom which is healthy, despite 
the appearance of higher derivative orders. Showing that the decoupling limit interactions of 
the helicity -0 mode are healthy despite the presence of higher derivative terms is therefore 
slightly more involved and the galileon dualities discussed in section 4 will allow us to do so 
without having to check order-by-order. 

3.3 Multi-Gravity 

Having covered the massive and bi-gravity cases, we now move on to multi-gravity models. 
For simplicity we will consider models which do have arbitrarily many spin-2 fields, but where 
all interaction terms only involve two such fields for any given term. Such setups have an 
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action 


S= d 


■D 


X 


3>i / D _ N 

Mpi ^ V~9{i) -^b(i)] + ^ | m \i,j)M pi \/~9(i) Y, P(ij),n Un{\J9~^ 9(j)) 


hi 


n =0 


(3.24) 


where we have once again assumed all Planck masses to be the same. Note that we have 
allowed for different coupling constants 'mup , although we could absorb this into the defi¬ 
nition of the coefficients which are now also allowed to vary from interaction term 

to interaction term. The corresponding decoupling limit action (as before after performing 
the Stiickelberg trick and projecting out the helicity-0 and -2 components) can be written 
as a superposition of linearised Einstein-Hilbert terms (EHL) and interaction terms ‘linking’ 
different spin-2 fields (for more on the link field interpretation see [28-32]) as follows 15 


?(*) _ 


__i M^cq/3 lM 

EHL ~ rl a^ 

A? 


odd) _ 

a link ~ 


3,0b) 


A3 




Stotai - y a i^s^hl+ y y 


j>i 


j(m) 
*link ’ 


h3 


(3.26) 


where (i,j) labels refer to the Stiickelberg fields resulting from the interaction term linking 
< 7 (i) and g^y As a result the dual to a field labelled by (i, j) is labelled by (j,i) - for details 
see section 4.1. Also note that each link has its own, independent set of coefficients fiy.j) 
implicit in the functions X and X. 


A concrete trimetric theory. As a concrete example consider the following trimetric 
theory 

S = j d D x ^My 2 ^g^R[g (l) ] + M°~ 2 ^^R[ 9{2) } + M°~ 2 R[g {3) \ (3.27) 

D _ D 

+ m M pi y/—g{ i) y^ Pn U n ^ yj <7^ <7(2)) + m M pi ^~9( 2 ) y^ cz n U n ^y/<7(2)17(3)) 

n=0 n=0 

Here we have three spin-2 fields with Einstein-Hilbert kinetic terms each. There are two 
bigravity like interaction terms with identical coupling constant m 2 , which makes taking the 

15 Technically our prescription for the decoupling limit now becomes slightly more involved. The limit we 
are taking here is 

m (i,j) —> 0, Mp —> 00 , Afixed, (3.25) 

where A 3 ^ p p Mpi. Note that we could have taken a different decoupling limit, where the coupling 

constants muj) scale differently and only some of the A 3 r i: y are kept fixed, while the others tend to zero. The 
essential point being that there is no unique decoupling limit isolating least suppressed interactions coming 
from each ‘link’ in the presence of several coupling constants in multi-gravity. 
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decoupling limit trivial, and we have two sets of constant coefficients j3 n and a n for the 
two links respectively (in the above notation these would be P n n 2 ) and P n .{ 2 ,S)i but we avoid 
unnecessary labels as much as possible for the concrete example here.) This theory is depicted 
in figure 1. The A 3 decoupling limit action now takes on the form 


S 






( 2 ) 

4 u (2) C ^ ,L a/3 


I uV-v caPfo(3) 
4 n (3) C ^ U a/3 


+ yC l (U*'“'(n ( i, 2) ) + 

+ Y h $ xl " { n (2l 3)) + • 


(3.28) 


We now demix this action at linear order via the linearised conformal transformations 


~> hfy + C( 1)2 )7r( li2 )?7, 

+ c ( 2 , l ) cr ( 2 , l)V + c (2,3) 7r (2,3) r ?) 

+ c (3,2) cr (3,2) r 7> (3.29) 

where the coj) are fixed by requiring the tensor-scalar interactions to vanish at linear order. 
This now leaves us with an action for the pure scalar piece of the schematic form 

S ~ J d D f 7r (l,2)b r n(n( 12 )) + 0-( 21 )C/ n (E( 21 )) + 7r(2,3)^n(S( 2)1 )) 

+ °'(2,l)b r n(n( 2 3 )) + 7T( 2 3 )C/ n (II( 2 3)) + <7( 3l2 )i/n(£(3, 2 )) I , (3.30) 


where we have suppressed constant (but n-dependent) coefficients for each term and the first 
and second line derive from the second and third lines in (3.28) respectively. 


Equations of motion. Let us again work out the eoms, this time for the trimetric action 
(3.30). We write 7 tv 1)2 ) = n and 7 i 7 2j3 \ = cf) and express the dual fields, a and p respectively, 
in terms of n and (j) (the same expression (3.18), and analogous for p in terms of (j>, still 
holds). From the bigravity case considered above we already know that irU n (H) + oU n {Y>) 
and (j>U n ($>) + pU n (P) give rise to second order Galileon eoms up to quintic order when 
expressed in terms of ir and (j). The interesting new piece comes from the scalar mixing 
generically induced by a multi-gravity theory, i.e. from the terms 

4>U n (Yj) and aU n (&). (3.31) 

~ ( 2 ) 

It is important to keep in mind that these terms come from ftp,} (Yip 2 ,i)) and h)xv -X' / “'(II/ 2)3 )), 
so we will give expressions in terms of the associated parameters (3 n and a n . The contribution 
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of (pU n ( S) to the 7 r eom up to quartic order in the fields in the action is 


£(1,7r) = 0 ) 

£(2,7r) = ^/5l < / ) a) 

£(3,77) = h {\(\>ab^ a + ^a^b) + @2 {^^ab - </>“aA), 

£(4,7r) = ^l(30afe C cVr a 7T 6 + \(t) a C b - \(f> a , a TT bc TT bc + 

+ @2 (4 > abc'K a 'K bC - (t> ab ^a C T^bc ~ ^ac^T^b + ^a^bc^ + ^^ab^c 

+ /3 3 (%4> ab lT a C ir b c - \4> a a TTb C T^ bC - \^> ab ^ab^ C c + f (t> a a^ b b^ C c) 

and to the (j) eom we find the following contribution 

£(2,0) = 

£(3,0) = -/3l(|7Ta67r a6 + 7;7r a 7rA) + /3 2 (i7T ab 7r a6 - 57T a a 7T 6 fe ), 

£(4,0) = MlvTaV^TTfee + TT a TT bc TT abc + ±7r a 7T a 6 7T 6 C c + iVr a 7T 6 7r a6 C c ) 

+ /3 2 (7r a a 7r fec 7T 6c - TT a C TT ah TT b c ~ ^TT^^abc + ^T^b^ac) 

+ /3 3 (|7r a c 7r a6 7r fec - f 7r a a 7r 6c 7r bc + |7r a a 7r fc fe vr c c ) (3.33) 

From cubic order onwards these are clearly higher-order in derivatives. In this sense the 
situation is worse than for the bigravity case, where even a disadvantageous choice of demixing 
procedure only resulted in higher-order derivatives from quartic order in the eom onwards. 
While we cannot jump to conclusions about the presence of a ghost-like dof here for the same 
reason as before (the Lagrangian is degenerate due to the non-local way in which 7r and <f> 
are related to a and p respectively), it is far from obvious that there are sufficient constraints 
in place to prevent extra dof’s from propagating. We will need the help of multi-galileon 
dualities, derived in the following sections, in order to show that this is nevertheless the case. 

The contribution of aU n {Q) to the 7r eom is (again up to quartic order in the fields in 
the action) 

£(l, 7 r) 2 dg, 

£(2,7r) = 2 d 0 ^ a a + § d 1 0 “ a , 

£(3,77) = d 0 (7r“ a 7r\ - TTabTT ab ) + a^^ab^ + § ^a^b) + 6t 2 {\(j> a atf b ~ ^ab ^), 

c _ ~ ,/2_ cab^ a _ be , 1 a b e \ 

^(4,7r) ^ ^be a^bc^ 3 ^ a ^ b^ c) 

+ dt 1 (l(j) a b c c ir a ir b + \(j) a c c^ a ^ b b - |^ a aVr fec 7T fec + \(t> a a^ b b^ c c) 

+ a 2 ((j) b b (l) a C c TT a - (j) bc (t> a b c lT a - \(t>ab(t> ab ^ C c + ^“a</> fc 6 Vr C c) 

+ ® 3 (g 4>a C (f> ab (frbe ~ \4 >a a4 > bc ( l )bC + b^ c) i (3.34) 
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The contribution to the eom, on the other hand, is 

c 3 * ^.CL 

£( 2 ,</>) 2 ^ 1 ^ a ’ 

^(3,0) = -di(|vr afe 7r afe + §7T a 7rA) + d 2 (</> a a A - (j) ab TT a b), 

£(4,tp) = d 1 (|7T a c 7r afc 7r fec + 3iT a 7r bc n abc + §7r a 7r a b 7r fe C e + § ir a ir b ir ab c c ) 

+ d 2 (0 a6 7r a c 7T fec - 0 a a 7T fec 7r 6c + (j) bc iT a ir abc - 4> b b ir a TT a c c ) 

+ Ol 3 (±(j) a C (j) ab 7T bc - \i> a a^ hC ^bc - \4>ab(l) ab ^ C c + l^a^b^c)- (3-35) 

Once again we have higher-derivatives dependencies from cubic order onwards and we will 
have to wait for the nrulti-galileon dualities, derived in the following sections, in order to show 
that these interactions are nevertheless healthy. 

In summary we see that, unlike for the massive and bigravity cases, the helicity-0 decou¬ 
pling limit of multi-gravity is not manifestly described by galileon interactions and thus also 
not manifestly ghost-free. The general proof of [5, 12, 14] for such models of course ensures 
that this limit must be ghost-free, but it would immensely help an investigation of the physics 
of such models to be able to write down the helicity-0 sector in a manifestly second-order 
form. This is the main aim of the remainder of this paper. In particular we have seen that for 
multi-gravity models the crucial difference to the bigravity case is the presence of additional 
cross-terms between different types of ‘dual’ and ‘original’ fields. This is a direct consequence 
of the ‘scalar mixing’ phenomenon discussed in [32]. In the following sections we will find 
that an extension of the single field galileon dualities [13, 22] will allow us to explicitly show 
that multi-gravity decoupling limit interactions are still of second order form and hence do 
not propagate any unwanted ghost-like dof’s, linking the multi-gravity decoupling limit to a 
particular set of multi-galileon-related interactions. 


4 Galileon Dualities 


We have already encountered the Sttickelberg scalar n and its ‘dual’ a in the previous sections, 
which are related by the non-local field redefinition (3.18). Here we briefly recap the origin of 
this duality in the bigravity picture [13], show that the duality is particularly straightforward 
to see in the ‘link field’ picture of of bigravity and collect some useful field relations, before 
extending the duality in the following section to the case of the multi-gravity theories (as 
introduced in 3.3). 


4.1 The Bigravity perspective 

Consider a generic bigravity theory 


S = d 


■D r 


i) ^[g(i)l + V~9(2) R[g( 2 )] + rn 2 V(g (1) ,g {2) ) 


(4.1) 
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where V is some potential interaction built out of the spin -2 fields g(i), g (2 ) and/or their 
inverses. Without the potential term this action is invariant under two copies of gen¬ 
eral co-ordinate invariance (two diffeomorphism symmetries), which we shall call GC^y and 
GC {2) . The potential term breaks these symmetries down to their diagonal subgroup. The 
Stiickelberg trick then transforms this action to one with additional fields (the Stiickelberg 
fields) and more symmetries (the full set of symmetries GC (p and GC ( 2 ) is restored). This 
symmetry restoration can be understood in the ‘link held’ formulation - for full details in the 
multi-gravity context we refer to [28-32]. 

In the bigravity context there are two possibilities to restore the full diffeomorphism 
symmetry of the action: Either, via the addition of (gauge) Stiickelberg fields, we transform 
g(i) into an object transforming under GC ( 2 ), or we transform g( 2 ) into an object transforming 
under GC^y This will make all the fields in the potential V transform under the same 
symmetry group - the potential will no longer break any symmetries, and as a consequence 
the full symmetry is restored. In terms of the Stiickelberg fields the two ‘options’ are 



or, in terms of functional composition notation, 

g(i) -> G( 2 ) = g(i) o y ( i >2 ), or g (2) -> G (1) = g (2) O y (2|1) . (4.3) 

In terms of the action (we only show the potential term, since the Ricci self-interactions are 
gauge-invariant and hence are invariant under the Stiickelberg replacement) this means the 
two possibilities correspond to 

Si = J d D x^/g^f( g(i),g( 2 ))-> J d D XyJg w oY { y 2 yg (2) ,) 

= J d D x^G^f(G (2 yg {2) ), (4.4) 

Sn = J d D x^/g^f(g w ,g^ 2) ) -» J d D x y /^ ) f(g [ 1 ) ,g { 2 ) oY { 2 y ) ) 

= J d D x^g^f(g {1) ,G {1) ), (4.5) 

Since the overall action is gauge invariant, we may perform a gauge transformation with 
parameter Y// 2 ) 011 or act with YY ^ on 5)” 4 , showing that 

W, = %,!)• 

y,i) = 1(1,2) (4.6) 

and we can now extract Stiickelberg scalars describing the helicity-0 mode for each ‘link held’ 
Y or Y (where we have implicitly set the helicity-1 modes to zero, as discussed above) 

y (l,2) -> *(1 )+d^(l,2), 

^(2,1) X (2) + ^ 1(T (2,1)i ( 4 - 7 ) 
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where we have suggestively named the second scalar field a. There are consequently two dual 
and physically equivalent ways of writing down the fully gauge invariant version of a bigravity 
theory such as (4.1), with the dual scalar fields n r and a. Note that, throughout this paper, 
we follow the sign convention of [22] for the relation between 7r and <7, which is opposite to 
that of [23]. 16 We emphasise that from the link field perspective the existence of this duality 
is a simple and direct consequence of the fact that there is no unique way to restore gauge 
invariance in the theory in question. Note that we have kept the labelling of the scalar fields 
as inherited from the link fields Y and Y = T -1 , so that the dual of 7T( 1j2 ) is < 7 ( 2 , 1 ), i.e. the 
label indices are reversed. 

We may now use (4.6) and (4.7) to implicitly define the dual field <7 in terms of 7r as 

(x + da) = (x + dir)- 1 . (4.8) 


We have already encountered this definition in (3.17) - the link formulation now makes it 
obvious where this definition comes from. One can now solve (4.8) in order to get an explicit 
expression for a in terms of 7 r. This relation can be written succinctly as Y(Y~ 1 (x))^ = x^. 
Expanding Y and T” 1 in terms of 7r and a respectively we then have 

(■ x + ditY + —— d u(x + dir) = x 11 , (4.9) 

d[x + dir ) M 

OO 

=Y d^a = -d^ir - ^D^d^a, where =—ir Vl ■ ■ ■ ir Un d ul ... Vn . (4.10) 

n =1 n - 

This final form allows us to recursively solve for d^a in terms of d^ir as 1 ' 

d^a =-d^(ir+~^ir u ir l/ + ...). (4.11) 

In (3.18) we already gave this expansion to quintic order. After some algebra one then reaches 
an expression for the n-th order piece in terms of lower orders (n is an order-label here, not 
a space-time index or a field label): 


n—1 


cr \_ n = - 


X D( n -i) [ a | 


2—1 


which can be solved to give 

OO - ^ 

+ 0<n. - IV X^ 1 ) 


n— 2 


a = —7T 


^ 2(n — 1)! 

n=2 v ' i=0 


n — 2 
i 


D {i) , 


(4.12) 


(4.13) 


where D^ n \(X) = d Vl ... Vn (n l/1 ■ ■ ■ ir Vn X) for some Lorentz scalar X. 18 We have now explicitly 
related the dual field a to 7 r. We again emphasise that in the link held picture the existence 

1(, These two conventions are related by sending cr —> —a. 

17 Note that the first terms in (4.11) are such that d^a is a total derivative and so it makes sense to talk of 
a. It is true [13], but certainly not obvious (since [D,d\ ^ 0), that this continues to all orders. 

18 We can immediately notice that (4.13) consists of Galileon terms and total derivatives, so cr appearing on 
its own in the Lagrangian is perfectly healthy when expressed in terms of n. 
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of the duality is a direct consequence of gauge invariance and hence is the result of a purely 
linear phenomenon (in a functional composition sense). 


4.2 Explicit field relations 

We can abstract the duality transformation away from the link field/multi-gravity argument 
above and simply view it as a field-dependent diffeonrorphism that treats a particular field, 
e.g. 7 r, in a privileged way. Here we collect some useful properties of this transformation from 
[23] and work out its form for the case of some particular ‘matter fields’, which will turn out 
to be relevant in section 5. 

We write the duality transformation linked to a field 7r as TV and now summarise its effect 
on all the objects in question. 19 As seen in the bigravity picture above, TV is a diffeomorphism 
at heart, so applying TV leads to a co-ordinate change 

TV : x » —^x^ = x^ + -^d'V(x), (4.14) 

where A is some scale (unsurprisingly we will find that A = A 3 in the bi- and multi-gravity con¬ 
texts discussed above), which we choose to set to unity in what follows. The Stiickelberg scalar 
7 T and its derivatives inherit their transformation properties from those of the Stiickelberg 
fields Y, so that in particular the derivative of 7 r transforms as a scalar under the action of 
TV (since Y ~ < 97 r) 


7r(x) ~ 

-A a(x) = 

dfj, 7 r(x) - 

-A d^a(x) 

, n-(x) - 

££(x) 


tt(x) - 4 (<9t r(x)) 2 , 

-<Vr(x), 

-[l“ + n“(x)]- 1 n2(a;), 


(4.15) 


where E^(x) = a{x) 20 and we have kept the dependence on co-ordinates x and x explicit, 

e.g. n is evaluated at x, £ at x here. Note that the transformation properties of n can be 
derived directly by taking the derivative of the second line in (4.15). When expressed in 
terms of the same co-ordinates (this amounts to iteratively solving (4.15)) then a and 7 r 
are related non-locally via (3.18), as discussed. The fact that we can also concisely express 
their relation in terms of (4.15) shows that this non-locality is a direct consequence of the 
(diffeomorphism) co-ordinate transformation mimicked by the Stiickelberg fields. In terms of 
n (or E) the Jacobian of the co-ordinate transformation x —> x induced by TV can now be 

19 We introduce the field-dependent label with an eye on the multi-field setups we will deal with in section 
5, where there will be a family of field-dependent dualities for each field. Note that this notation is somewhat 
different from that used in [23] . [23] are dealing with a single field galileon duality family, so use of a field label 
would be redundant, but they do use a further parameter s labelling different types of duality transformations 
in this setup. We have just taken s = 1 here for simplicity, absorbing any s into the definition of A, but an 
explicit generalisation to arbitrary s is straightforward. 

20 In index-free notation we may also write the expression for E as E = — [1 + II] -1 II or equivalently as 
E" 1 = -(l + ir 1 )- 
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succinctly expressed as 


= det(l + II(x)) 1 = det (1 + £(x)) , (4-16) 

assuming that the sign of the determinant is positive . 21 As far as fields other than n are 
concerned, is just a diffeomorphism, so a scalar %, vector V or arbitrary-index tensor T 
transform just like under diffeomorphisms 


dx a 

~dxb 


: < 


X(x) x{x) = x(x ), 

y^x) -> Vn(x) = j£iVu(x) = [1 + E(®)]£ Vu(x ), 
T^ n {x) T^ n (x) = [1 + E(®)]£ ■■•[! + E(®)]£ T Ul .' Mri 


(4.17) 


(®)> 


Note that we have chosen to write the hnal expressions in terms of £(x), but we could have 
just as well written them in terms of n(x), e.g. Vfi(x) = [(1 + ’\\.{x))~ 1 \ v (] y v {x). The associated 
inverse duality transformation V~ l simply involves swapping a and n in the above mappings 
- for details we refer to [23]. Here it will turn out to be very useful to work out the explicit 
action of V n on another scalar x and also on d^x and d^d^x- We then have 


X(x) ->• x{x) = x(x ), 



dnX(x) -> duxix) = || ?d„x(x) = [1 + 'Z{x)] v lt d IJ ,x{x ), 
d^d v x{x) d^x(x) = [1 + E(x)]“ d a ([1 + E(x)]£ d/9x(®)) ■ 


(4.18) 


Throughout the remainder of this paper we will frequently use these maps and talk about 
‘applying’ duality maps or duality transformations to an action or eoms. It is worth empha¬ 
sising that by this we mean e.g. expressing V v (x) in terms of V v {x) by solving the mapping 
(4.17) for V v (x) in terms of V v (x). It does not mean simply replacing V v (x) —> V u (x). 


4.3 (Single) Galileon dualities 

A general n-th order Galileon interaction term for a single held <r may be written 

SgV) = J d D xcrU n (Ti) . (4.19) 

If we view a as the dual held for 7 r, we may now explicitly substitute for a via (4.13) and, 
after some algebra, hnd 

J d D xoU [n) {Y,) = \yi) n {n+l) J d D x f kl (p D n _ n fc + iy ^ U (n+k-i)(ir) + d^jA , 

(4.20) 

21 For a discussion of this requirement see [22]. 
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which we see is precisely of the form ‘Galileons + total derivatives’ (the exact form of J^ is 
not important here), and upon partial integration is equivalent to 


j d° x a U{n] (Y.) = (-1)" +1 („+1 )fd D x (jt {k + 1)( c _ (“)!(* - ^ ' 

(4.21) 

A single field Galileon in terms of a can therefore be expressed as a Galileon in terms of n 
by direct substitution. In particular this means that a scalar field action as arising in the 
decoupling limit of Hassan-Rosen bigravity 


S 


D 


d° x C( n ) (tt C/( n) (n) + crf7 (n ) (£)) 

n=0 

D 

d D xY b (k) 7TU (k)( U ) 


k =0 


(4.22) 


manifestly describes the do/of a single held Galileon (and is hence ghost-free) 22 , where the 
coefficients bi n ) are given by 


fr(n) 


n 

c («i + 


k=0 


(—l) fe+1 (A; + 1)(D - k)\ 

(n + 1)(D — n)\(n — k)l ^ 


(4.23) 


Instead of explicitly substituting for a and subsequently rearranging the action into 
Galileons plus total derivatives, we may equivalently follow [22, 23] and instead use the ab¬ 
stracted duality transformation maps discussed above in 5.1 to map a Galileon in one frame 
into its dual form. Under the duality transformation T> n we then have 


5 Gal (7T) = f d D xY^U {n) (U) (4.24) 

** n 

- J dD x ^2 c (n) det(l + E(x)) ^o-(x) + ^5 M o-(x)a^cr(x)^ U {n) 

which after renaming the dummy integration variable x into x and using some of the algebraic 
properties of the characteristic polynomial U/ n ) can be brought into the form 


~£(x) 

1 + £(x) 


cGal/ \ P-rr 
S (n) W 


D 


d D x V C(„, (- l)" +1 (A; _ ( ° )!( / ) i fc)! (^ + £(*)], (4.25) 


k=n 


22 Note that because Y(Y 1 (x)) = Y 1 (Y(a;)), the bigravity picture neatly ensures that these relations take 
exactly the same form with n and a swapped. 


- 23 - 









where all variables are functions of x now. Now we can already pick out crU^ (£) as a Galileon 
term. For a^a^U^(E) we may integrate by parts to find 


d D xa^U {n) {E) = (D - n)\ / d D xa^ 


■ ay 


2 (D-n)\ 
n + 2 
2 (D - n) 
n + 2 




d D xcT[/ (n+1) (S) oc 5g^ 1} (<r) 


Gal 


Upon substituting this into (4.25), up to boundary terms we arrive at 



(- 1 ) 


n+1 


D 

(n + 1) / ^ 

/c=n 


(D-n)! 


(fc + l)(D- jfe)!(A: 




(4.26) 


(4.27) 


which manifestly describes a Galileon in terms of a. As expected this has precisely the same 
form as (4.21). 


5 Multi-Galileon Dualities 


Knowing how to map single held Galileon interaction terms from one held picture into its 
dual form, we now move on to multi-Galileon interaction terms. A general such interaction 
term for N helds at order n in those helds may be written [41-43] 


eM 

‘'multi—Gal 


N 


M2,-,A 


[M2,-,Mn] 7r (H) 7F fe) 


M2 

M 2 


i\ —1 


7T, 


(in) 


Hn 


(5.1) 


where the (ij) are indices labelling different scalar helds and aw 1 ’"'>M are constant coefficients. 
For the purposes of the remainder of this paper we will restrict ourselves to investigating 
dualities relevant in the context of the multi-gravity theories discussed in section 3.3 above. 
As a result two-held Galileons of the following type will be of particular relevance for us 


cW 

°bi—Gall 


= “ (1 ' 2 . 2| C 


Mn 

•,Mn 


1^(1 ■ 


vr, 


( 2)1 


( jD _ n + l)i 7r(1)C/ («- 1) ( n (2)) • ( 5 - 2 ) 


The most general multi-gravity theories naturally will include more general multi-Galileon 
type interactions in their decoupling limit - we leave the explicit proof of this to future work. 


5.1 Demixing and the multi-gravity decoupling limit 

Multi-gravity theories as discussed above in section 3.3, i.e. theories with arbitrarily many 
spin-2 helds and bigravity-like interaction terms, in the decoupling limit have the following 
interaction terms (for helicity-2 and -0 modes) 


5 = -- 


E + 5 E (AU+J+^n^)) + Al i( y ) ++'“'(%i,)) . (5.3) 
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In order to demix this action at linear order, one can employ the following transformations 


hfy -> h® + Y ( c (*J) 7 r (*,i) + Viu,, ( 5 - 4 ) 

j 

where the constant coefficients and duj) are fixed by the requirement to eliminate 

mixing between the helicity-2 and -0 modes at linear order, they are in general a function 
of the particular coupling constants (m(t -y /3 n ,(i,j)) chosen for each bigravity-like interaction 
term in the original action, and we insist on coj) = 0, if d( t jj ^ 0 and vice versa 23 and 
C(jj) = d(ij) = 0 if there is no interaction linking gu\ and - for details see [32], As a 
result, in the demixed decoupling limit, we end up with the pure scalar piece of the action 
taking the following form 


j^i k^i,k^j / \ 

+ yi yi ( c {j,k)'K(j,k)U n (^‘(j,i)) + d(i t k) a (i,k)U n (J[(ij)) 

i,j k \ J 

j^i k^i,k^j / \ 

“k ^ ] y ] ^3,{i,j) ( c (*,fc) 7r (i,fe)^n(!!(*,j)) T 

i,j k \ / 


(5.5) 


The first line describes Galileon-like interactions for 7Toj) and the dual fields just as we 
encountered in the bigravity case. The second and third lines describe new cross-interactions 
between sets of different 7r’s and cr’s that are qualitatively different. Also note that we have 
no interaction terms of the type cr(i t k)U n {^(k q) and - these are eliminated 

by choosing a demixing procedure like M\ in (3.19) that yields non-zero 
and interactions, i.e. the first line in (5.5). The third line describes cross¬ 

interactions, which are already explicitly Bi-Galileons in terms of 7r’s or cr’s. Finally, it is 
worth pointing out that, for some simple theories, this structure simplifies further, e.g. for 
line theories (see figure 1) one can make a consistent choice of and d^ ^ that makes the 
third line in (5.5) vanish. 

In that action we have N fields nuj) as well as their duals V(j.i) ■ However, since this 
action describes only N fundamental scalar dof’s, we would like to express the full action 
(and assess features such as its ghost-freedom, interaction scales etc.) purely in either the 7r- 
or a- frame. 21 Here we choose to map the full action into the a -frame (this choice is of course 
completely arbitrary). From section 4.1 above we know how to map the self-interaction terms 
from one frame into another and that a single field Galileon in one frame maps into a single 
field Galileon in the other. However, here we now also have to deal with the new cross-terms, 

2J This choice is always consistent and corresponds to the demixing field re-definition for each h( t \ being a 
function of either iru j) or ou j) for all j, but never a function of both. 

24 Technically we could choose one field from each 7r — n-pair - which field is labelled 7r and which a for each 
link is purely conventional. 
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whose mapping we will investigate here. In total we would like to establish the following 
mappings 


7r (1) tf (n) (n (1) ) ^ 

a (2) u (n) { n (1) ) ^ 

vr (1) t/ (n) (S (2) ) ^ 
7r ( 2)C/( re )(n (1) ) ^ 
7r (1) C/ ( „)(n (2) ) ^ 


A)(C(1)), 

Cr (2) ) 5 

CT(2))i 

~ 1)^(2) 

^///(o‘( 1) ! 7!'(2)) -> ^///( cr (l)> cr (2)), 

~ X^7T ^2) 

^/v(°'(l)) 7r (2)) -•4/V'( <T (l)> cr (2))> 


(5.6) 


figuring out explicit expressions for all the A in the process and in particular confirming that 
the dual formulation of a given term describes precisely the same number of dof’s as the term 
being mapped. In other words we would like to explicitly check that the dynamics of the 
fields is still described by healthy eoms and consequently that it is explicitly free of ghost-like 
instabilities. From section 4.1 we already know that 

* / ^o( (7 (l)) 'y ' J c (n)V(l)U(n)(^-‘(l))i (5-7) 

n 

so we will proceed to investigate the cross terms now. Note that ^( 2 )U( n \ (II(i)) and (11(2)) 

are of course already explicitly Galileons in the 7r-frame, but we are here interested in how 
terms map to the cr-frame. 


5.2 Dualities I: Mapping the action 

We begin by investigating the mapping of cross-terms in (5.6) using duality maps at the level 
of the action, in direct analogy to section 4.3. 


Bi-Galileon dualities I: The first term we will consider here is the following cross¬ 
interaction term 


®"(2) ^(n) (H(i)) 



1), CT(2)). 


(5.8) 


We now aim to employ duality transformations to fully map this into er-space at the level of 
the action, i.e. to find *4/(cx(i),ov 2 )). Comparison with (5.2) shows that (5.8) is a two-field 
Galileon in terms of the two fields <T( 2 ) and 7T(i). Here we will show that this is also true when 
expressed in terms of a/ 2 ) and C(i). Taking a general superposition of terms like (5.8) and 
applying the duality transformation Thqi) we find 


S(5V( 2 ),7r (1) ) 


J d D xa {2) U {n) (U {1) ) 



d D x det(l + S (1) (x))<T (2 )I7( n ) 


-S ( i)(x) 

1 + S(i)(x)_ 


(5.9) 
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where it is important to note that <j( 2 ) transforms like any other ‘matter’ scalar field according 
to (4.18) under T>^^y As far as a^ 2 ) is concerned we are applying a diffeomorphism to the 
action and, it being a scalar, <j( 2 \ remains invariant under that transformation. After renaming 
the dummy integration variables x to x and some algebra as before, (5.9) can then be brought 
into the form 




5 §) 1 ( <7 (2)>’r (1) ) ^ (-l) n / d 


D, 


D 

T. 

k=n 




(5.10) 


which is manifestly a two-field Galileon in terms of dm and <7( 2 ). In other words, for some 
constant coefficients C( n ) defined via (5.10), we have 


1),C(2)) = y ^C(n)(T(2)U(n) [ S (l)]~ 
n 


(5.11) 


Bi-Galileon dualities II: Moving on to the second cross-term in (5.6) we now consider 


7r (l)t / '(n)(S(2)) 



•A.h{(T( 1),U(2)). 


(5.12) 


We could straightforwardly map this term into 7r-space using the above recipe, but here 
we want to uniformly map all possible terms into c-space, so we again apply the duality 
transformations and find 


S(n)W),CT (2) ) = 



X 



x 'y ] c (n) 7I "(l) (S( 2) ) 

n 

i D x^2 c (n) det(l + S (1) )(d (1) + 

n 

i( 1 + E (i)r 1 ]°9o([( l + s (i)r 1 ]>f)! 


(5.13) 


We may now further manipulate this action, e.g. we can extract the factor of [(1 + E^))^ 1 ]” 
from the elementary symmetric polynomial U( n ) in order to cancel the determinant factor 
in front, but the essential feature here is already visible at this stage: This is an action 
explicitly higher than second order in derivatives. The dependence on means we have 

a dependence on third derivatives of an) and we may confirm that the eoms obtained from 
varying (5.13) with respect to and <7( 2 ) are indeed higher order in derivatives as a result 
- we will do so below. 

We consequently have a non-Galileon interaction term here and the corresponding eoms 
will be higher order as well. However, this does not mean that we necessarily have an Ostro- 
gradsky ghost, since as discussed above in section 2.2 a system can still be secretly second 
order even though it possesses higher order eoms due to the presence of constraints. This is 
most easily seen at the level of eoms , so below we will move on to investigate the multi-Galileon 
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dualities at that level. As expected, this will turn out to be much more straightforward for 
the second cross-term considered here than mapping at the level of the action. 


Bi-Galileon dualities III: Before doing so, let us briefly comment on the last two 
terms in (5.6). As far as 'Dtt(i) is concerned 7r(2) and transform in exactly the same way, 
i.e. as a scalar field under diffeomorphisms, so the final two terms in (5.6) transform precisely 
like the previous two. As such we have (from (5.10)) 

[ d D xY,c {n) n { 2 ) U {n) [U {1) } V -^ f d D x ^ ^T(—ir C(n) 7T( 2) I/ (fc) [S (1) (x-)], 

J n J n k=n ^ M >’ 

(5.14) 

and hence 


1),7T(2)) = y^C M TT (2 )U(n) P(l)]- 


(5.15) 


This can now be fully mapped into the <r-frame via (5.13), losing its manifest Galileon form 
in the process as before. For the final cross-term in (5.6) we find (just as in (5.13)) 


d ° X yi c (n) 7r ( l) u (n) t n (2)] 


V 


”■( 1 ) 


x U, 


& D x y2 c {n) det(l + S (1) )(<T (1) + 


(n) 


[(! + S (l)) ( [C 1 + S (1)) i ]^ 7r ^ ) 


(5.16) 


Again this can be fully mapped into the c-frame using Pw 2 p which will not, however, restore 
manifest Galileon form. An(a^,a^ 2 )), i), 072)) an< 4 Aiv(.&( i), cr ( 2 )) are therefore not 

of Galileon form and we will now move on to working at the level of the eoms in order to show 
that, despite their higher-derivative appearance, these dual actions describe healthy theories. 


5.3 Dualities II: Mapping the equations of motion 

We start by considering a generic multi-galileon, whose action is given in (5.1), which we will 
here schematically write as 

5 (n) 1 ( 7r (i)> 7r (2)i---i 7r (n)) ( 5 - 17 ) 

The corresponding i eoms that arise from varying this action with respect to each ir^ are 


Sr 


(i) 


II II II 

(1) ’ (2)’ • • ' ,7F W 


= 0 , 


(5.18) 


where Roman numerals label the maximal order in derivatives. So, for example, n 1 denotes 
that the eom depends on first derivatives of it. In the particular case considered here, multi- 
galileon eoms, the eoms of course only depend on the via their second-derivatives as is 
required by the Galilean symmetry (i.e there is no direct dependence on the fields themselves 
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or on first derivatives of the fields). 


Bi-Galileon duality dictionary: We now go back to our dictionary of Bi-Galileon 
terms (5.6), for which we wanted to establish their mapping under duality transformations 
and schematically work out their eoms and how these eoms transform under the duality. Our 
aim is to show that, when completely mapped to the a- frame, the resulting eoms can be 
put into an explicitly second-order form, demonstrating that no additional ghost-like dof’s 
propagate. If, for a given set of terms, we can show that purely second-order eoms in one 
frame transform into purely second order eoms in the other frame, then we have successfully 
shown that Galileons in one field representation map to another set of manifest Galileon 
interactions under duality transformations - as we saw in the previous section this does not 
have to be true for all types of multi-galileons. 25 We have 


7T(i)[/ (ri )(n(i)) = 

^ £ 

^ C 7T(1) 

CT (2)t / (u)(n ( i ) ) = 

^ £ 

^ C 7T(l) 

7T (l)U(n){'£(2)) = 

=b. £ 

^ C ^(2) 

7r (2)t / (n)( n (l)) = 

^ C 7T(1) 

7r (l) t/ (n)(n(2)) = 

^ £ 

^ °^(2) 


.// 

(i) 


II 

(1) 

II ' 
a (2) 

and £ CT(2) 

'ii ' 

_ 7r (f)_ 

II 

(1) 

II' 

°(2) 

and £ n(J) 

'll 

( 2 ). 

II 

(1) 

II 

7r (2)_ 

and £ 7T(2) 

. 7r W. 

II 

(1) 

7F (2). 

and £ 7T(]) 

7r (2) 


(5.19) 

(5.20) 

(5.21) 

(5.22) 

(5.23) 


In other words, and again as required for multi-Galileon interactions, the eoms only depend 
on 7T(j) and via II(j) and respectively. Let us now go through term by term and es¬ 
tablish the transformations of these eoms when duality transformations are applied in order 
to map all fields into the a- frame. 


Galileon dualities: First the single-field Galileon interactions (5.19) 


^(1) [^(l)] £tt {1) [ (1 + ^(1)) ^(i)] — & 


'( 1 ) 


ii 

a (i) 


(5.24) 


so that the duality transformation manifestly turns a purely second-order equation for 
into a purely second-order equation for This proof is of course equivalent to the one via 
mapping the action shown in section 5.2, which proved that a single Galileon action in terms 
of 7r straightforwardly maps into one in terms of 


Bi-Galileon dualities I: Now for the first cross-term (5.20) with eoms 


£, 


( 1 ) 


II 

^(1)^(2) 


and £, 


( 2 ) 


II 

^(1) 


(5.25) 


25 In other words, multi-galileon eoms in one frame can map into seemingly liigher-order eoms in another 
frame, but still be physically equivalent due to the presence of constraints. 
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Remembering the transformation properties of IR) and S( 2 ) ((4-15) and (4.18) respectively), 
under these eoms are mapped to 


£, 


a) 


ii 

^( 1 )^( 2 ) 


V 


7T(1) 




( 2 ) 


7r 


1 T> 


7T(1) 






(i) 


( 2 ) 


II ^ ~ IJ 

a (l)i a (l) > Cr (2)> Cr (2) 


a 


(5.26) 


These eoms are higher-derivative in 07 ^ now. 
differentiate to obtain 

£ 

C <T( 2 ) 


iii 

°{i) > 


ii 

a (i) 


However, we may take the second eom and 




(5.27) 


We can now use £ a{0) to solve for a^ in terms of cr^ and insert the solution back into (5.26). 
In that way we finally obtain a set of eoms which are explicitly second-order 


£, 


■(i) 


TI I II 
^(l)) (T ( 2 )’ CT ( 2 ) 


and £ c 


( 2 ) 


a 


(5.28) 


This demonstrates that, in the a-frame, a cross-term in the action like (5.20) does not lead 
to an Ostrogradsky-ghost via higher-derivatives in the eoms (i.e. signalling a dependence on 
additional initial conditions). Note, however, that in working with the above schematic form 
for the eoms, we have not used all of the information available. In particular, while (5.28) 
shows that the eoms were secretly second order and hence free of any Ostrogradsky-ghost, it 
does not show that they are purely second-order, i.e. manifestly of the (multi-)Galileon type. 
We know that this is the case from (5.10), i.e. am and 07 2 ) in fact obey Galilean invariant 
eoms, so they are indeed bona-fide Bi-Galileons in the cr-frame. However, without explicitly 
solving (5.27) for and carrying out the explicit substitution, the eom analysis carried out 
here would not show that the dependence on a ^ in (5.28) vanishes. 


Bi-Galileon dualities II: We now look at the eoms coming from the second set of 
cross-terms (5.21) 


£r 


( 1 ) 


a 


and £ c 


( 2 ) 


II 

7T(1)^(2) 


(5.29) 


We once again perform the duality transformation T > 1T(1) and find 


Sr 


( 1 ) 


II 


V 


7T(1) 


S e 


'( 2 ) 


II 

^( 1 )^( 2 ) 


1 V 


7T(1) 


Sr 


Sc. 


( 1 ) 


( 2 ) 


II ^11 

a (l)’ CT (l) ’ a (2) 


ii ill r ii 

a (l)’ a (l) > Cr ( 2 )’ Cr ( 2 ) 


(5.30) 


This time we have a dependence on in both eoms. Now, we could find the linear 
combination of these eoms that eliminates the third time-derivative of am (and hence, by 
Lorentz-invariance, also a 1 ^ as a whole), solve the resulting equation for <j(i) (the solution 
will be given in terms of two initial conditions for am and two for and then solve one 

of the original equations for cj( 2 ) (since we have already solved for the fact that there 

is a third-derivative dependence on 071 ) in the final equation does not introduce extra initial 
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conditions). However, there is no obvious way of eliminating the higher-derivative dependence 
in (5.30) and obtaining two equivalent eoms which are explicitly second-order in analogy to 
(5.28) (at least in the absence of specifying further information about the form of the eoms 
(5.30)). 

Nevertheless, there is an alternative way of showing this. We may take (5.29), use the 
second eom to solve for cr( 2 ) and insert that solution into the first eom. As a result we will 
obtain a set of different, yet physically equivalent, differential equations schematically of the 
form 


ii 



and £-2 


7 I fi)’ 


ii 

CT (2) 


(5.31) 


Expressing the dynamical evolution of the system via the equations £\, £2 we can now proceed 
just as for (5.25) above. £\, £2 are of precisely the same form as (5.25), showing that our 
system secretly obeys second-order eoms in the cr-frame of the form 


£ 


A) 


11 1 11 

C7 (1)J °(2p °(2) 


and £ c 


( 2 ) 


II 

^( 1 ) 


(5.32) 


This explicitly shows that terms in the action like (5.21) do lead to secretly second-order 
eoms and hence do not lead to Ostrogradsky ghosts in the cr-frame. Note, however, that in 
this dual formulation of the explicitly Galilean-invariant eoms (5.29), the Galilean symmetry 
is not manifest just as in the previous example. 


Bi-Galileon dualities III: To fully map the remaining terms into the cr-frame we need 
to apply both and P 7r ( 2 )- Still at the level of the action, using (5.10) we can map 


d X7T ( 2 ) C/ (n) (n (1) ) - > 2^ d ® £ (n) 7r (2) C/ '(n)( S (l))> 

n ** 

^7T (2) % ^ I 

d X 7 T ( 1 ) [/ (n) (n (2) ) - > 2^ d x£ (n) 7r (l) C/ (n)(S( 2 )), 


(5.33) 

(5.34) 


where the C( n ) are constant coefficients as defined via (5.10). The eoms from (5.33) are 
consequently 

and £ a crfL irfl ■ (5.35) 




( 2 ) 


II 

^(l) 


We now apply V n ( 2 ) to these eoms and, in precise analogy to the mapping of (5.29) to the cr- 
frame, find that in the cr-frame all solutions are captured by the following set of second-order 
differential equations 


Si 


11 1 11 

<T (2)> f7 (l)) CJ (l) 


and £ 2 


The eoms from (5.34), on the other hand, are 


£r. 


( 2 ) 


II II 


and £ T] 


(i) 
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CT ( 2 ) 


II 

a ( 2 ) 


(5.36) 


(5.37) 


Applying T> n ^ to (5.34) and proceeding as above, in the cr-frame all solutions are captured 
by 


£\ 


„ii „ 1 „ii 

a {l)-> a (2)' a {2) 


and £2 


a 


(5.38) 
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Again we have shown, at the level of the eoms, that interactions in the action (5.22) and (5.23) 
do not lead to Ostrogradsky ghosts in the a-frame. However, the associated eoms , which were 
explicitly of Galileon form in the 7 r-frame, do not display manifest Galilean symmetry in the 
a-frame (again, just as in the previous example). 

5.4 Dualities III: An explicit multi-galileon example 

The above discussion showed that, while we can use the duality mappings to prove that multi¬ 
gravity theories (3.24) are free of any (Ostrogradsky) ghosts in the decoupling limit when 
mapped into either the a- or 7r-frame, the resulting interactions and eoms are not always 
manifestly of Galileon form. While the duality is a direct consequence of diffeomorphism 
invariance (as the link field picture shows) and the fact that the o’/'/r-relation is an invertible 
field re-definition guarantees that the same dof’s are described in either frame, it can in 
principle be difficult to read off symmetries of the helicity -0 decoupling limit interactions. 

For this reason we now briefly give an example of a multi-metric theory in which one 
can explicitly show that only multi-Galileons appear in the decoupling limit and the second- 
order nature and a Galilean shift symmetry are therefore manifest, i.e. we do not have 
the complications of (secretly second-order, but) seemingly higher-order eoms. Recall the 
trimetric theory (3.27) and the corresponding demixed decoupling limit action (3.30), which 
we will here schematically write as (suppressing n-dependent constant coefficients) 

<5 ~ / d D X^2 ( 7T(1,2 )Un [11(1,2)] +°'(3,2)^n [^(3,2)] 
d n 

+ ( cr (2,l) + 7r (2,3)) fin [^(2,1)] + U n [n (2 ,3)]) ^ , (5.39) 

where the first two terms are inherited from mixing with and /i( 3 ) in the pre-demixed 
action (3.28) respectively (they come from the nodes at the ends of the graph), and terms on 
the second line are inherited from mixing with fy 2 ) (tluey come from the middle node). 

From (5.7) we see that the first two terms, and two arising form the final one are manifestly 
of Galileon form and hence ghost-free when written either in terms of 7 r* or cr*, leaving the 
two cross-terms terms: 


tt( 2 , 3 )U n ( s ( 2 ,i)) and (J( 2 ,i)U n (n (2 , 3) ) . (5.40) 

Simply using (4.13) in either yields a series of multi-galileon terms as well as a series of 
additional terms the lowest order of which would be d a TTid a TTi£™(irj), which is a Galileon for 
i = j, but leads to higher order equations of motion otherwise. Now it is if course possible 
to use the arguments of sections 5.2 and 5.3 to demonstrate the healthy nature of the whole 
theory, but in this case it actually possible to approach it in a different way and explicitly 
show that only multi-Galileons are present. 

Looking at the cross terms we immediately note that they only involve <7( 2 ,i) and vt( 2 , 3 ) , 
and not tt(i, 2 ) or a ( 3 , 2)1 and when written in terms of these fields are manifest bi-Galileons and 
hence healthy. Thus if we declare that our two ‘fundamental’ helicity-0 degrees of freedom 
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Figure 2. Different orientations of the link fields (as indicated by the arrows, and labels above each 
edge) lead to taking different fields as the ‘natural’ helicity-0 dof’s (indicated below each edge). In order 
to obtain manifest multi-Galileon interactions for hclicity-0 dof’s in the trimetric theory considered 
here, the key if for links to be directed away from only nodes of degree one. 

are 07 2 ,i) and 7r ( 2 , 3 )? not n (i, 2 ) and n ( 2 , 3)1 these terms are now manifestly healthy bi-Galileons. 
The remaining terms will always be healthy single-field Galileon interactions regardless of 
which fields are considered ‘fundamental’, where it is worth re-iterating that for each pair of 
7 tuj) and <?(j,i) • which field is considered to be the ‘fundamental’ one and which the ‘dual’ is 
purely a matter of convention. Using only single-field dualities we may therefore schematically 
write (5.39) as 

s ~ J d D X ^2 ^ a ( 2 ,l)U n [^( 2 ,l)] + 7T(2 j 3)C/n[n(2 j 3)] 

+ ^"(2,3) Uji [S( 2j i)] + (J(2,l)Un [11(2,3)]^ j (5-41) 

which as an explicit bi-Galileon interaction for the two scalar helicity-0 dof’s (T( 2 ,i) and tt( 2 , 3 ) - 
In terms of the multi-gravity theory this has a nice interpretation in terms of the direction 
of the (Stiickelberg) link fields mapping a metric form one site to another. The original 
theory has three metrics <7(i),(2),(3)? and expressing all interactions in the 7r-frame amounts to 
introducing link fields Yn 2 ) mapping g ( 2 ) onto site (1) in the first interaction term (connecting 
g (1 j and g ( 2 )) and Y( 2 , 3 ) mapping g ( 3 ) onto site (2) in the second interaction term (connecting 
g( 2 ) and g ( 3 )), see figure 2. Expressing all interactions in the a -frame would amount to 
reversing both links. If, however, we solely change the direction of Yn 2 ) so that it instead maps 
5 ( 1 ) onto site (2), then we are naturally poised to express things in terms of {<T( 2 ,i) 5 tt( 2 , 3 )}> 
again see figure 2. 

This approach works nicely in the trimetric case considered here, and will continue to 
work in a theory for which there is only one node of degree greater than one, such as a star 
graph [20, 32]. However in other cases it fails, due to the fact that terms such as <JiU n {Yij) are 
not manifestly of ghost-free Galileon form when written in terms of 7Tj and ttj. In such cases 
we must rely on the arguments of sections 5.2 and 5.3 to demonstrate the healthy second-order 
nature of the interactions. 
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6 Adding a coupling to matter 

Throughout this paper we have dealt with Multi-Gravity theories of the type 


S = 2 / dPxy/-g^R [s (i) ]+m 2 M£j 2 / d D x^/-5 (1) F (5 {1) ,.. ■, 5(JV)))+5matter[^, 

i=l ^ •* 

( 6 . 1 ) 

where the matter action <S m atter[$i> {<7(j)}], describing the way all matter dof’s 4>j couple to 
gravity, has been ignored so far. 

Minimal coupling and dualities: The standard, minimal and GR-like, coupling to 
gravity in Massive, Bi- and Multi-Gravity theories amounts to minimally coupling all matter 
dof’s to a single metric only 

^matter j = ^matter > (6-2) 

which we choose to label with the label index (1) here. If the matter metric only couples to one 
other metric (carrying a label (2)), i.e. the matter metric is a node of degree one in the theory 
graph, then the matter action will only be sensitive to a single helicity-0 mode it. Expanding 
around a flat background space-time r]^ v and focusing on the interactions between that single 
helicity-0 dof and the matter dof’s, from now on for simplicity taken to be described by a 
single scalar dof \ (no important feature depends on this assumption, however), we have 

^matter = J d D xC ma ,tte r \x( x )i 9 nX( x )> K(x), 8^ 1 T(x), 8^8"TT(x)\ , (6.3) 

where we have assumed that there is no dependence on second derivatives or higher of the 
matter fields (which is straightforward to generalise via (4.18), however). Under the duality 
transformation this maps to 

x{x),[{l + T 1 )-\d v x{x), (6.4) 

- (°( x ) + \{dcr{x)) 2 ^j , -d^a(x), -[(1 + T l )~ 1 } v lx d v d a a(x) , 

where we have replaced the dummy variable x with x (i.e. after performing the duality trans¬ 
formation). Unlike for single-field Galileons the form of matter interactions is therefore not 
invariant under the duality (as was also not the case for some of the multi-Galileon interac¬ 
tions considered above). However, both ways of writing down the interaction are of course 
physically equivalent due to the nature of the duality transformation, and the associated eoms 
will consequently remain secretly second-order in the sense discussed in the previous sections 
(also see [23]). 

Demixing and Decoupling I: We stay with the case where the matter metric only 
couples to one other metric (carrying a label (2)), as is the case in standard Bi- and Massive 


qi) 


^matter 


d D x 


det(l T S(.T))/2 ma tter 
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Gravity. As part of the demixing procedure we will then perform (up to choices between 
different consistent demixing procedures, e.g. (3.19)) 

+ C(l, 2 ) 7 r (l, 2 )^- (6-5) 

As a result we find that the decoupling limit matter action is 

Smarter,dec = J d° XTr(x)T^ \x(x), dx(x)] , (6.6) 

where Tjf is the trace of the stress-energy tensor defined in the conventional way with respect 
to metric g^y This is straightforward to see, if we choose to introduce Stiickelberg fields 
solely through the metric(s) to which matter does not couple, since the only dependence of 
the matter coupling on helicity-0 modes then comes in through (6.5). Under the duality 
transformation this decoupling limit contribution maps to 

Smarter,1,dec = J d°XTr(x)T^ [x(x), dx(x)] (6.7) 

- J d D x det(l + X(x)) ^cr(x) + ^(<9cr(x)) 2 ^ T" [y(x), [(1 + S) _ 1 ]^y(x)] , 

where we have replaced dummy variables as before and we note that, even though the decou¬ 
pling limit contribution only depends on n and not its derivatives, in the dual picture there 
is explicit dependence on a, da, d 2 a. 


Demixing and Decoupling II: If the matter metric couples directly to M other metrics 
(i.e. it is a node of degree M in the theory graph - see figure 1), then expanding around a 
flat background space-time 77 ^ and focusing on the interactions between that single helicity -0 
dof and the matter dof’s as before, we find the following schematic action 


^matter 


= / d D xC 


matter 


xO), d/lX(x ), 7 T( 1 ) 0 ), 0 * 0 , (x ), 

• ■ •, 7 T(M) 0*0 , dfjTT (M) (x ), d^d u ir {M) (x) 


( 6 . 8 ) 


The demixing procedure for the matter coupling node (1) now takes on the following form 
(as above we still assume matter only couples to a single metric) 

M+l 

+ ^2 ( c (i,j) 7r (i,j) + d(i,j) a (i,j)) *7/mo ( 6 -9) 


3 =2 


leading to a decoupling limit action 

. M 

Smatter,dec ~ I d X ( 7r (n)(®))^[x 0*0, <9*0*0] 


( 6 . 10 ) 


n =1 


V 


7T(1) 


j d D x det(l + S ( 1 ) (x)) |- ^cr ( 1 ) (x) + ^(da^x)) 2 ^ + ^ (vr (n) (x)) 

xTZ[ X (x),[(l + X m )-XdvX(x)], 


- 35 - 




where we have used a shorthand notation for the helicity-0 fields 7 rn j) = n(j-i) to avoid 
clutter and ignored constant coefficients throughout. In the same fashion one can now itera¬ 
tively apply the remaining T> v ^\ to fully map this into the cr-frame. Once again the duality 
guarantees that even if the resulting eoms are higher order, the initial value problem is well- 
defined in terms of two conditions for each scalar field and hence no extra (Ostrogradsky) 
dof’s propagate. 


Non-standard matter couplings: Until now we have only considered matter couplings, 
in which matter minimally couples to a single metric. However, in general we may have a 
matter coupling of the form 


^matter[*1*1 > — ^matterj Sfiu [l7(l) j ■ • • j ff(JV )]]j (6.11) 

i.e. we couple to more than one metric, where we have insisted that the weak equivalence 
principle is upheld and matter minimally couples to an effective matter metric g™ ttcv that 
is a function of the other metrics. This can be done consistently (in the decoupling limit) 
following the matter couplings proposed by [51, 52], As a result the matter action, when 
expanding around a flat background space-time ri tw and focusing on the interactions between 
that single helicity-0 dof and the matter dof’s as before, takes on the form 


^matter 


= / dPxC 


matter 


x(x), dftX(x),7T W (x), dft r (1) (x), 


• ■ •, 7T(JV) (x ), (x ), d^d u 7 r (J v) (x) 


( 6 . 12 ) 


Investigating duality mappings in the decoupling limit for these new matter couplings is left for 
future work. However, one may reasonably expect that these mappings will be rather different 
from the cases considered above in one important respect. A primary reason why the matter 
coupling (6.2) was not invariant under the duality map was that (6.2) explicitly breaks the 
symmetry between metrics present at the level of the potential self-interactions. Introducing 
Sttickelberg fields via one metric or the other in a bigravity-like interaction, while physically 
equivalent, no longer leaves the form of interactions invariant at the level of the matter action. 
The many metric coupling (6.11), especially when considering maximally symmetric matter 
couplings in multi-gravity theories [52] , can restore this symmetry and hence potentially some 
of the duality invariance as we encountered it for single-held Galileons. We leave answering 
these questions for future work. 


7 Summary and Conclusions 

In this paper we derived the decoupling limit for a particular class of multi-gravity theories 
and investigated the resulting interactions for helicity-2 and helicity-0 modes. By making use 
of the link-held picture, we have shown how this is directly related to the existence of dual 
descriptions for such theories. In particular, we extended the known Galileon dualities to a 
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set of multi-Galileon dualities, relating a class of multi-Galileons to other multi-Galileons and 
additional higher-derivative interactions, which are healthy due to the existence of constraints. 
These dualities were discussed in several complementary ways: At the level of the action, 
the equations of motion, as related to diffeomorphism invariance and as abstracted field re¬ 
definitions/mappings. Finally we also showed how matter couplings transform under duality 
transformations depending on the nature of the coupling. 

Understanding the decoupling limit of multi-Gravity aids us in investigating the physics 
of these theories and their (low energy) interactions. An obvious task for the future will 
be to complete the work started here and derive the full decoupling limit interactions for 
multi-Gravity (i.e. for all ghost-free such models and for all helicity dof’s ) and to use this 
to understand the cosmological phenomenology and relevance of such theories. The dualities 
uncovered in the process are also of interest from a purely field-theoretic point of view, es¬ 
tablishing the equivalence of seemingly unrelated field theories via invertible, non-local field 
re-definitions that would have been very hard to discover in another way. Realising that 
these dualities exist can be extremely useful in investigating multi-scalar field models, e.g. 
for setups where the duality relates strongly and weakly coupled theories and hence makes the 
physics in seemingly strongly-coupled regimes calculable (this is also the case for the single 
Galileon duality). It also alerts us to the existence of a rich space of higher-derivative theories, 
which are nevertheless healthy and free of any ghost-like degrees of freedom. Investigating 
theories of interacting spin-2 fields appears to be an especially fruitful avenue to discovering 
healthy theories of this kind. 
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A Appendix: Quintic order Lagrangians and quartic order eoms 


In section 3 we computed the decoupling limit helicity-0 eom for ghost-free Bigravity up to 
quartic order in the fields in the Lagrangian and hence cubic order in the eom. Here we 
complete this task up to quintic order. Extending (3.22), where the eom was worked out after 
linear interactions had been demixed with Adi from (3.19), £( 5 ) is given by 


%) = /3 1 (-VT a C 7r a6 7r ; , d 7r cd + \^ a a^^^cd. + \'KabK ab 'Kcd'K Cd ~ 7T a a TT b b TT cd TT Cd + g7r“ a 7T b fe 7r C c 7r d d ) 

+ P2d' K a C TT ab 'K b d lTcd ~ 6n a a TT b d TT bC TT cd - 17T. ab TT ab 7T cd 7T Cd + §71° a TT b b TT cd TT Cd ~ |tt“ a 7T b b TT C c TT d d ) 

+ P 3 (-6Tr a c 7T ab TT b d 7r cd + 8TT a a TT b d 7T bc 7r cd + 3ir ab TT ab ir cd TT cd - 67r a a 7r b ;,7r cd 7r cd + n a a 7T b b 7T C c TT d d ) 

+ /3 4 (|vr a c 7r a6 7rfc d 7r cc ; - ™7r° a 7r 6 'V c 7r cd - |vr ab 7r ab 7r cd 7r cd + §7r a a 7r b b 7r C(i 7r cd - ^n a a TT b b n c c iT d d )■ 

(A.l) 
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Noticeably this is still completely second-order in derivatives, so manifestly free of Ostrograd- 
sky ghosts. Extending (3.23) on the other hand, where demixing had been performed with 
M .2 (3.19), we find 


£(5) = /3i(-7j7T a C 7r“ 7 T b 7 T cd + g7T a a 7 T b 7T C TT cd - j7T a a 7T b TT cd TT C + ^7T a a 7T b 7T C c TT d 

- lTT a TT b d TT bC TT acd - \'n a TX bc 'K bC 'K a d d + ±7r a 7r\7T C c 7r a ' i d - ^TTaTT^TThcd ~ ^T^TXa^TThcd 

- lTT a TT a b TT b C TT c d d - \-R a TT b 7T , ab TT/ d ~ ^“tT^TT^TT ab cd + jVbVVcTT ab d d - \l\ a TXa^hc d ) 

+ P 2 {\TXa C lX ab TX b d TX cd ~ 37t“ a TT b d TT bC TT cd - TT ab TT ab TT cd TX Cd + 27t“ a 7T b b TT cd TT Cd - \TT a a TT b b TT C c TT d d 

i o_a _ d^bc^. /Tr d __b _cc£_ rrr a m-b —-C _ d | _a_ b^cd^. , 1 ^a^b^ cd _ 

i 7T& 7T 7T ac d 7T 7T ^7T 7 T ac d 7T 7T fr7T c^a d ^ a 7 T 7 Tfocd ' 2^ ^ ^bed 

- TT a TT a b TT C c TT b d d ~ \ir a TT b TT a C c TT b d d + 7r“7T b 7r Cd 7T ab cd ~ ^^TT 0 c 7T af) d rf) 

+ P 3 (—^TX a C TX ab TX b d TT cd + 6n a a Tr b d TT bC TT c d + §7T ab TT ab 7T cd 7T Cd - §7r“ a 7r fe 6 7r C(i 7r cd + § TT a b TT C d ) 


+ /5 4 (|t r a C 7T a TTfc 7r cd - ^7T a a 7r 6 TT C TT cd - §7T ab 7T a 7T cd 7r“ + §7T a a 7T b 7T C(i 7r C - ^7T a a 7T 6 7r C c 7r d ). 

(A.2) 


As had already been the case for £/ 4 \ for this choice of interactions, there now is explicit 
higher-order dependence on derivatives, now inside both the /3i and /?2 terms, rather than 
just inside the /3i terms as was the case for Interactions are of course still healthy 

following the logic of section 2.2, but this serves to emphasise that the choice of demixing 
procedure can be helpful in making the true number of propagating dof’s manifest. 
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